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LIMIT  THEOREMS  FOR  PERIODIC  QUEUES 


1 . Introduction  and  Summary 

Let  A = (A(t),  t > o)  be  a non-stationary  Poisson  process  with  bounded 
measurable  arrival  rate  function  x(*)  satisfying  X(n  + t)  = x(t)  for  all  t 
in  [0,  1 ) and  nel  = (l,  2,  . . .}  .We  call  the  unit  of  time  a day,  and  we 
refer  to  a point  t in  [0,  1)  generically  as  a time  of  day.  Defining 


X = / X(t)  dt 


we  say  that  A is  a periodic  Poisson  process  with  a period  of  one  day  and  an 
average  arrival  rate  of  X .We  denote  by  T^,  T^,  . . . the  jump  points  of 
the  right -continuous  process  A 

Let  , Sg,  ...  be  independent  and  identically  distributed  (i.i.d)  posi- 
tive random  variables,  independent  of  A , with  distribution  function  (d.f)  F 
having  F(0)  = 0 . We  define  o = XE(S^)  and  assume  throughout  that  o < 1 
For  t > 0 let  X(t)  = + . . . + , with  the  convention  that  an  empty 

sum  is  zero,  and  then  let  Y(t)  = X(t)  - t and 


Z(t)  = Y(t)  - inf {Y(s-):0  < s < t}  . 


We  interpret  T^,  Tg,  ...  as  the  customer  arrival  times  at  an  ordinary  single 
server  queueing  system  and  S^,  Sg,  ...  as  the  customer  service  times.  Then 
Z = (Z(t),  t > o]  represents  the  server  load  process  (or  virtual  waiting  time 
process)  for  this  system,  assuming  that  the  server  is  initially  idle;  cf.  Benes 
(1958)»  Alternatively,  one  may  interpret  Z as  the  contents  process  for  a dam 
with  input  process  X = (x(t),  t > 0)  and  (constant)  unit  release  rate,  assum- 
ing the  dam  is  initially  empty.  If  the  queueing  context,  the  actual  waiting  time 
of  customer  n is  given  by  Wn  = Z(Tn~)  = Z(T^)  " sn  • is  t*ie  purpose  of 


( 

this  paper  to  study  the  asymptotic  behavior  of  the  processes  Z and  W = 

{W  , ne  I)  . 
n' 

Queues  with  non -homogeneous  Poisson  arrivals  have  been  studied  earlier  by 
Takacs ( 1955 ) , Reich  (1958,  1959)  anc*  Hasofer  (196k).  Hasofer  (196b)  considered 
the  case  of  periodic  input  and,  imposing  some  further  conditions  on  the  arrival 
rate  function  and  service  time  distribution,  showed  that  the  probability  of 
server  idleness  is  an  asymptotically  periodic  function  of  time. 

In  Section  2 we  focus  attention  on  the  integer  time  points  (i.e.,  where 
days  begin)  at  which  the  system  is  empty,  showing  that  these  times  constitute 
a sequence  of  regeneration  points  for  the  continuous  parameter  process  Z 
A similar  sequence  of  regeneration  points  is  identified  for  the  discrete  para- 
meter process  W . In  Section  3>  these  results  are  combined  with  discrete 
renewal  theory  to  prove  the  following  weak  limit  theorems.  For  each  time  of 
day  te  [0,  1)  there  exists  a proper  distribution  such  that 

P{z(n  + t)  < x)  -»  Ht(x)  as  t — > 00  for  all  x > 0 . Furthermore,  there 
exists  a proper  distribution  G such  that  pfw^  S x)  -*G(x)  as  n ->  00  for 
all  x > 0 . We  also  prove  strong  limit  theorems  for  Z and  W . Speci- 
fically, for  all  x > 0 , 

t 

d)  tf l{z(s)  < x}ds  "H(x)  as  t->°° 

0 

with  probability  one  (w.p.l)  and 


As  a by-product  of  (2)  we  shov  that 


1 

;(x)  = 


) ds  , x > 0 


Finally,  in  Section  ^ we  show  that  the  limit  distributions  G and  H are  re- 


lated by 


= 1 - p+  p*(G*F)(x)  , x>0  , 


where  * denotes  convolution  and  the  d.f.  F is  defined  by 


E{S1)  -F(x)  = /[  1-F(y)]dy  , x>0 


In  particular,  this  shows  that  H(0)  = 1 - p , which  is  more  or  less  obvious 
from  physical  considerations.  The  relationship  (3)  between  the  asymptotic  vir- 
tual and  actual  waiting  time  distributions  is  exactly  the  same  as  that  found  by 

Takacs  (1963)  for  the  Gl/G/l  queue.  Our  line  of  proof  follows  Lemoine  (197M 
/ 

and  shows  that  Takacs'  result  actually  holds  for  much  more  general  systems  if 
one  interprets  the  distributions  H and  G as  sample  path  limits  as  in  (l) 
and  ( 2 ) . 

The  results  presented  here  are  weak  in  the  sense  that  we  neither  compute 
nor  characterize  the  distributions  G and  H , except  to  show  that  each  can 
be  obtained  from  the  other  via  (3).  By  demonstrating  that  interesting  asymp- 
totic distributions  do  exist  for  periodic  systems,  however,  we  hope  to  suggest 
some  potentially  tractable  problems  in  the  theory  of  queues  with  non-stationary 
Poisson  input  and  hence  to  rekindle  interest  in  this  important  class  of  models. 


2.  The  Regenerative  Structure 


We  formally  assume  a probability  space  (n,  P ) on  which  is  defined  a 
right -continuous  compound  Poisson  process  X = (X  (t),  t > 0)  with  unit  junq> 


SVMUe; **vwrmw**em 


rate  and  jump  size  distribution  F.  Let 


A(t) 


t 


X(u)  du 


0 


t > 0 . 


Then,  a(*)  is  non-decreasing  and  continuous  with  a(o)  = 0 and  A(n  + t)  = 
nX  + X(t)  for  nel  and  t e [0,  1)  . Let  X(t)  = X*(A(t))  for  t > 0 , 

and  let  A(t)  be  the  number  of  jumps  of  the  process  X = (x(u),  u > 0}  during 
[0,  t]  . Thus,  A and  X are  described  in  Section  1.  Also,  let 


■9*  = &(  X*(u)  , 0 < u < t)  , t > 0 , 


0<t<<» 


and 


& = '^A(t)  * ® — * — 00  * 

We  say  that  a random  variable  T : fl  -»  fO,  »]  is  optional  if  (T  < t}  e for 
all  0 < t < oo  . And,  for  T optional  we  define 

^ = {E  : Efl  {T  < t}  e.Wt  for  all  0 < t < ®)  . 

The  process  X is  strong  Harkov  (cf.  Blumenthal  and  Getoor  (1968),  pp.  37-bh), 
from  which  we  easily  obtain  the  following. 

Proposition  1 . If  T is  optional  and  finite  w.p.l,  then 


jirewai^^fc.  .<***»* 


i 

t 

: 


I 


f 


I 


pjx(T+ t)  - X(T)  = 0|«^Tj  = expj-[A(T  + t)  - a(T)| 

w.p.l  for  t > 0 

Proposition  2 . If  N is  integer-valued,  optional,  and  finite  w.p.l,  then 
(X(N  + t)  - X(N),  t > 0}  is  independent  of  and  distributed  as  {X(t), 

t > 0)  . 

We  define  Y(t),  Z(t),  T and  W in  terms  of  X as  in  Section  1.. 

Let  O.  = inf  (ne  I:Z(n)  - 0}  , with  QL  = <»  if  the  indicated  set  is  empty. 
Since  Y(t)  and  Z(t)  are  ^-measurable  for  t > 0 , it  follows  that  a 
is  an  integer-valued  optional  time.  Recall  that  we  assume  p = \E{S^)  <1  , 

and  observe  that  E{X(l)}  = p . 

Proposition  3 • E {a } < a> 

.g- 

Proof:  Since  Y(t)  = X (A(t))  - t , we  have  Y(t)/t  -»p-l  as  t-»» 
(w.p.l),  and  hence  Y(t)  -*  - w as  t -*  » (w.p.l).  From  the  definition  of 
Z it  then  follows  that  sup{t:Z(t)  = 0)  = » w.p.l.  Moreover,  it  is  clear 
that  sup(t:Z(t)  > 0}  = oo  w.p.l.  Let  = 0 and  recursively  define 

V = inf (t  > | : Z(t)  > 0) 
n n ' 

and 

§ , = inf {t  > V : Zft)  = 0) 

n h l n 

for  nel  . In  queueing  terminology,  g is  the  time  at  which  the  nth  idle 
period  begins,  and  is  the  time  at  which  the  nth  busy  period  begins, 

nel  . From  the  remarks  above,  the  variables  t and  V are  finite  w.p.l 
for  nel  , and  each  is  clearly  optional.  Let  = inf  {k  e I : k > | ) for 
nel  , so  that  7^  is  the  first  integer  time  following  the  onset  of  the  nth 


■I 

5 

i 


j 

5 


5 


period.  Then  a = yx,  where  N = inffne  I : y < v ) . Now. 

n n — n 


P{N  > n + 1)  = P{N  > n,  v < 7 ) 

— — ’ n 'n 


= P {a  > ln,  x(yn)  - X(|n)  > 0} 


P{X(7n)  - X(|n)  > 0|®n)dP 


(a>  | ) 

— n 


where 


9 = JT  . Since  y < | + 1 , we  have 

n | n — bn 


p{x(rn)  - x(in)  > o|&  )<  p|x(fJn  + l)  - x(^)  > o|«  j 

= 1 " exp|- [A(ln+1)  - A(in)]| 


by  virtue  of  Proposition  1 . But  A(t  + 1)  -A(t)=X  for  all  t > 0 , and 

ao  we  have 

P(N  > n + 1)  < (l  - e"X)p(a  > |n)  - (l  - e"X)p{N  > n)  . 

n -1 

Thus,  by  induction,  P(N  > n)  < (l  - e ) for  nel  , which  implies  that 
E{n)  < eo  . in  particular,  a = 7N  and  so  we  have  a < °o  w.p.l  . Note  that 
the  total  idle  time  of  the  server  during  [0,  Ct]  is 

N-l 

- inf {y(s-)  : 0 < s < a)  = ^(vn  - $n)  + ( 7 „ - !„)  • 


Hence,  from  the  definition  of  Z , we  have 


0 


N-l 

= z(a)  = Y(a)  +]^vn  ~ + (rN  " - Y(a)  • 

n=l 


But.  the  nth  idle  period  has  duration  (v  - £ ) < 1 for  1 < n < N - 1 , 

and  7jj  “ Ijj  S ^ as  well>  80  that  0 < Y (a)  + N . Thus,  we  have 

0 < -Y (a)  < N w.p.l,  and  P(-Y(a)  > o}  ;>  P[A(1)  = 0)  > 0 , lienee 

0 < E{-Y(a)}  < E{N)  < « . 

Having  shown  p(qj  < °°}  = 1 , we  set  = a and  then  recursively  define 

0^+1  = inf  jn  e I : n > and  Z(n)  = 0 j 


for  kel  . It  follows  from  Proposition  2 that  Y(a^)  , Y(a2)  ■ Y(a^)  , 
. . . are  i.i.d.,  as  are  a^,  - a^,  . . . . We  write 


(5) 


-YK> 

k 


k 


Then,  letting  k-»oo  in  (5),  we  conclude  from  the  strong  law  of  large  num- 
bers that  E{-Y(a)3  = (1  - p)E{aJ  . Thus,  (1  - p)E(a)  < E[n)  < <»  , which 

completes  the  proof. 


Remark:  If  p > 1 , one  can  show  E(al  = « by  observing  that  a is  a 

weak  descending  ladder  index  for  the  random  walk  (Y(n),  nel}  . 

Throughout  the  remainder  of  the  paper,  we  define  the  integer-valued 
optional  times  {a  , nel}  as  in  the  proof  of  Proposition  3,  and  we  set 
= A(&n)  for  nel  with  & = 0^  . Since  a is  optional  and 
(A(n)  - A(n  - 1),  nel}  is  an  i.i.d.  sequence,  Wald’s  Lemma  gives  E(&)  ~ 
E Cot } • E (A ( 1 ) ] = \E(a]  < co  . Also,  setting  = (3  s 0 , we  define 


°Tl 

Vx)  = fhz{t)<  x}dt 

an-l 


and 


< x) 


for  n e I and  x > 0 , with  the  convention  that  an  empty  sum  is  zero. 

(Thus,  vn(x)  = 0 wittl  positive  probability  for  each  fixed  n and  x.)  Let 

U(x)  = U^(x)  and  V(x)  = V^(x)  , and  note  that  U(x)  < a and  V(x)  < p 

Since  an  is  optional,  it  follows  from  Propositions  2 and  3 that  (Z^^-t-u)  , 

u > 0}  is  independent  of  Jf  and  distributed  as  Z . The  following  is 

n 

then  immediate. 

Proposition  h . Fix  x > 0 . The  sequences  {an  - an_j,  nel)  , 

(p  - p , , nel}  , (U  (x),  nel}  and  {V  (x),  nel)  are  each  i.i.d. 
with  finite  mean. 


3.  The  Limit  Theorems 

For  each  te  [0,  1)  we  define  a proper  d.f.  Ht  by 


E(a)  ' Ht(x) 


a-1 


k=0 


' {z(k  + t ) < 


x > 0 


Proposition  5 • For  each  t e [o,  1)  and  each  x > 0 , 


1 

4 


< 


’i 

i 


•x 


i 


I 


I 


3 


8 


lim  P(z(n  + t)  < x}  = Hfc(x) 


Proof:  We  have  already  observed  in  Section  2 that  {Z(a  + u),  u > 0}  is 
independent  of  & and  distributed  as  Z . Thus  for  n - 0,  1,  2,  . . . 


p(Z(n  + t)  < x)  = Pl'a  > n,  Z(n  + t)  < x) 


11 

^ P{a  = k}  • p[z(n  - k + t)  < x) 


for  fixed  t e (0 , 1)  and  x > 0 . Let  = p{z(n  + t)  < x}  , 
bn  = P(a  > n,  z(n  + t)  < x}  , and  = p{a  = n}  for  n = 0,  1,  2, 
Then,  the  equation  above  can  be  equivalently  expressed  as 


(6)  vn  « bn  +]£fkVn-k  , ® - 0,  1,  2,  . . . . 

k=l 

Now,  since  f^  > 0 and  E fa)  < <®  it  fo>.  ows  from  the  discrete  renewal 
theorem  (cf.  Feller  (1968),  p.  550)  that 


00 

lim  Vn  =]Cbk/E 

n"*“  k=0  / 


1{a  > k)1(z(k  + t)  < x)|/E^a^ 


the  second  equality  holding  by  virtue  of  monotone  convergence.  And,  the 
the  second  expression  for  the  limit  is  clearly  Ht(x)  , completing  the  proof, 


H(x) 


w.p.l 


Proposition  6 . For  all  x > 0 


(z(s)  < x] 


ds 


as  t ->  oo  where 


H(x) 


0 


Proof:  From  Propositions  3 and  If.  and  the  strong  law  of  large  numbers  it 

is  easy  to  show  that 


1 

t 


t 


(Z(s)  < x} 


ds 


e {u(x)  )/e  (a)  w.p.l 


I Hi  liiTrtigaiV^-*  "* 


the  last  equality  following  fromFubini's  Theorem  (cf.  Neveu  (1965),  p.  91), 
Comparing  this  with  the  definition  of  H^(x)  completes  the  proof. 


Proposition  7 • For  all  x > 0 


lira  P{W  < x]  = G(x) 

n — v ' 

n — i<» 


where 


G (x)  = E (V(x) )/  E O ) 


Proof : Since  a is  optional,  we  see  that  {w^  + n>  nel3  is  independent 

of  .?  and  distributed  as  W . Thus 


P(W  < x)  = Via  > n,  W < x]  +y^  P(B  = k3  P {W  , < x) 

n — — n — / j n - k — 


for  fixed  x > 0 and  nel  . Let  v = p{w  < x)  , b = p{ft  > n.  W <x}  , 

— n n—  7 n — 7 n — 7 

and  f = P{f3  = n)  for  n = 0,  1,  2 . . . where  V = b_  = 0 . Then, 

n 7 7 7 7 0 0 

the  equation  above  has  the  form  (6)  with  f^  > 0 , and  again  we  invoke  the 

discrete  renewal  theorem  to  obtain 


CO 

lim  P{W  < x)  = V bR  / E { 
n — >00  , ^ / 


e Ev.^v.x)  H 


11 


h 


1 * .. 


I , - , ~ • -f 'f  .wm.-r-Nv-.- 


The  last  expression  for  the  limit  clearly  equals  E(V(x))/E{pj  , completing 
the  proof. 

Combining  Proposition  U with  the  strong  law  of  large  numbers  and  the 
definition  of  G , we  easily  obtain  the  following  result. 

Proposition  8 . For  all  x > 0 


k=l 


(W,  < x] 

k — 


G(x)  w.p.l 


as  n — > oo 


1 

Proposition  9 . X G(x)  = J " x(s)  Hg  (x)  ds  , x > 0 


Proof:  Fix  x > 0 and  let 


r(  t)  = 


A(t) 


k=l 


(W,  < x] 
k — 


for  t > 0 . From  Proposition  8 and  the  fact  that  A(t)/t  -»  X w.p.l  as 

t -» oo  , it  follows  readily  that  y(  t)/t  -»\G(x)  w.p.l  as  t ->  » . Further- 
more, 0 < y( t)/t  £ A(t)/t  , and  E(A(t)/t)  = A(t)/t  < 2X  for  all  t > 1 

Thus,  uniform  integrability  gives  E(7(t)/t}  -*  \G(x)  4s  t ■*  » . Now  let 


! 


« 

i» 

i 

i* 


: 31 


► 

i4 


n-l 


rn(t)  = Xrf  1(z(kt/n)  < x)  [A^k  + 1)t/“)  -A(kt/n)J 
k=0 

for  nel  and  t >0  . Then  7 (t)  ->  /(t)  w.p.l  as  n — » 00  , and  hence,  by 


dominated  convergence,  Efy^t)}  — > E Ct'C t ) j as  n 


From  the  independent 


1? 


increments  of  X we  have 


} 


x > 0 


\ 

\ 


U (x) 


a 

1 Cz (t ) > x) 


dt 


Then  U(x)  + U (x)  = a for  all  x > 0 , where  U(x)  is  defined  as  in 
Section  2.  Furthermore,  U(0)  and  U (0)  represent,  respectively,  the  total 
server  idle  time  and  the  total  server  busy  time  during  the  initial  regenera- 
tive cycle  [0, a)  • Thus  we  have 


A (a) 

(7)  a = u*(0)  + 0(0)  = ^ sk  + u(o)  , 

k=l 

where  Sk  is  the  height  of  the  kth  jump  of  X (the  service  time  of  the 
kth  arriving  customer)  and  we  follow  the  convention  that  an  empty  sum  is 
zero.  Now,  from  (3)  in  Lemoine  (197*0  it  then  follows  that 

A(a) 

(8)  U(x)  = ^min|sk,  (x  - Wfc)  j + u(0)  . 

k=l 

Using  (7)  and  (8),  we  thus  have 

A(a)  A(a) 

- Q - »(->)  - E sk  - L Yk  ■ 

k=l  k=l 

WR)+|  . Using  Wald’s  Lemma  we  can  obtain 


E(A(a)}  ■ E { s x ) = XE(a)  •E(s1)  = pEfct)  . 


lit 


1 

i 


Let  denote  the  o -algebra  in  .?  generated  by  T , and  then  for 

kel  let  ctf,  denote  the  cr-algebra  generated  by  T , S,  , . . . , T S 

1 1'  ’ k*  k* 


and  T,  . . Then 

k + 1 


E EK 1 


I k A{A(ct)  > k - 1}! 


- Ee!eKk-i)  • 1 


{A(cO  > k - 1} 


!A(a) 

EE(YkK-l) 

fc=l 


j A(a) 

E ) 2 E{si)  *E[(x  " V+] 


Proceeding  exactly  as  in  the  argument  for  (5)  in  Lemoine  (I97U),  it  follows 


( A(a) 

(u>  * E\ 

( k=l 


= E(S1)  -E{A(a)}  * (G*F)(x)  = pEfo)  ’(G  *F)(x) 


Since  E{u(x)}  - E{a)  * H(x)  , the  desired  result  follows  from  (9)>  ( 10 ) and 
(11).  This  completes  the  proof. 


mean  oi  F is  E{S^  }/?E(S^}  , so  that  Proposition  10  gives 


*■*•*?■  mm******':-*-*?  ? -'■ 


* • ^ * • '--'■  ■*v 


x H(dx) 


X E {S  ^ } / 2 + 


5 . Concluding  Remarks 

In  defining  our  model  (Section  2)  we  have  assumed  A( • ) to  be  non- 
decreasing and  absolutely  continuous  with  a periodic  density  x(*)  , but 

the  assumption  of  absolute  continuity  is  not  really  necessary.  We  may 
instead  assume  we  are  given  a non-decreasing  and  continuous  function  A 
with  a(0)  = 0 , a(°0  =»  , and  A(n  + t)  = nX  + A(t)  for  nel  and 
te  [0,  1)  . All  of  our  results  continue  to  hold  exactly  as  stated,  except 

that  \(s)ds  is  replaced  by  A(ds)  in  Proposition  9. 

This  suggests  the  further  generalization  where  A is  permitted  to  have 
discontinuities.  In  that  case,  the  customer  arrival  process  A may  have  jumps 
of  any  integer  size  at  time  points  where  A is  discontinuous.  (The  size  of 
the  jump  in  A at  such  r.  time  point  is  Poisson  distributed  with  mean  equal  to 
the  height  of  the  jump  discontinuity.)  Our  results  for  the  server  load  pro- 
cess Z continue  to  hold,  but  those  for  the  waiting  time  process  W must 
be  altered  somewhat . 

With  A continuous,  our  model  can  also  be  generalized  by  allowing  the 
service  time  distribution  to  depend  on  the  time  of  day  at  which  a customer 
arrives.  With  the  traffic  intensity  p properly  redefined,  all  of  our 
results  except  Proposition  10  continue  to  hold,  but  the  proof  of  Proposition 
3 becomes  much  more  complicated. 
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